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■ Abstract 

cd ■ 

c/3 , In the equilibrium statistical mechanics of classical Coulomb fluids, 

the long-range tail of the Coulomb potential gives rise to the Stillinger- 
Lovett sum rules for the charge correlation functions. For the jellium 
model of mobile particles of charge q immersed in a neutralizing back- 
I ground, the fixing of one of the g-charges induces a screening cloud of 

Q ' the charge density whose zeroth and second moments are determined 

O . just by the Stillinger-Lovett sum rules. In this paper, we generalize 

these sum rules to the screening cloud induced around a pointlike guest 
charge Zq immersed in the bulk interior of the 2D jellium with the cou- 
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pling constant T = Pq {/3 is the inverse temperature), in the whole 
region of the thermodynamic stability of the guest charge Z > —2/T. 
The derivation is based on a mapping technique of the 2D jellium at the 
coupling r = (even positive integer) onto a discrete ID anticommuting- 
[ field theory; we assume that the final results remain valid for all real 

I values of T corresponding to the fluid regime. The generalized sum 

■ rules reproduce for arbitrary coupling T the standard Z = 1 and the 

trivial Z = results. They are also checked in the Debye-Hiickel limit 
\ r — > and at the free-fermion point F = 2. The generalized second- 

^ • moment sum rule provides some exact information about possible sign 



oscillations of the induced charge density in space. 
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1 Introduction 



The present paper deals with the equihbrium statistical mechanics of a classi- 
cal (i.e. non-quantum) jellium, sometimes called the one-component plasma, 
formulated in two spatial dimensions (2D). 

The jellium model consists of mobile pointlike particles j = 1, . . . , N of 
charge q and position vectors r^, confined to a continuous domain A. The par- 
ticles are embedded in a spatially uniform neutralizing background of charge 
density —qn. The bulk regime of interest corresponds to the limits N ^ oo 
and |A| ^ oo with the fixed particle density n = A^/|A|. 

According to the laws of 2D electrostatics, the particles can be thought 
of as infinitely long charged lines in the 3D space which are perpendicular 
to the confining 2D surface A. Thus, the electrostatic potential at a point 
r G A, induced by a unit charge at the origin 0, is given by the 2D Poisson 
equation 

A(/)(r) = -27r(5(r). (1.1) 

For an infinite plane A = R^, the solution of this equation, subject to the 
boundary condition V0(r) as |r| — > oo, reads 

0(r) = -ln(^), (1.2) 

where r = |r| and the free length constant tq will be set for simplicity to 
unity. In the 2D Fourier space defined by 

fir) = /^/»exp(ik-r), (1.3) 
f d^r 

f{k) = J — /(r)exp(-ik-r) 
the Coulomb potential (11.21) exhibits the form 

m = Y2 (1-5) 

with the characteristic singularity at A; = 0. This maintain many generic 
properties of "real" 3D Coulomb fluids with the interaction potential 0(r) = 
1/r, r e R^. 

Because of the presence of the rigid background, the equihbrium statistics 
of the jellium is usually studied in the canonical ensemble under the condition 
of the overall charge neutrality. The 2D statistics depends on the coupling 
constant F = (3q^ with (3 = l/{k-BT) being the inverse temperature; the par- 
ticle density n only scales appropriately the distance. Let the symbol (■ ■ ■)/3 
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denotes the canonical averaging. At the one-particle level, one introduces the 
average number density of particles 



n{r) = {J2S{r-r,)) . (1.6) 

At the two-particle level, one introduces the two-body density 

n(2)(r,r') = (E'^(r-r,)5(r'-rfe)\ . (1.7) 

It is also useful to consider the pair correlation function 

hir^r') = , , ' - 1, 1.8 
^ ' ^ n{r)n{r') ' ^ ^ 

which tends to at asymptotically large distances |r — r'| — > oo. 

The bulk jellium is in a fluid state for high enough temperatures, i.e. the 
density of particles is homogeneous, n{r) = n, and the two-body density is 
translation invariant, n^^)(r,r') = n*^^^(|r — r'|). There are indications from 
numerical simulations [Ij that around F ~ 142 the fluid system undergoes a 
phase transition to a 2D Wigner crystal. In what follows, we shall restrict 
ourselves to the fluid region of F-values. 

Through a simple scaling argument, the exact equation of state for the 
pressure P, j3P = n[l — (r/4)], has been known for long time [2]. The jellium 
is completely solvable, like any Coulomb system, in the high-temperature 
Debye-Hiickel (DH) limit T — >^ [3], characterized by a monotonic exponential 
decay of the pair correlation function h{r) at asymptotically large distances 
r ^ oo. The systematic F-expansion of statistical quantities around the DH 
limit can be done within a bond-renormalized Mayer diagrammatic expansion 
The 2D jellium is mappable onto a system of free fermions at the special 
coupling r = 2 [5]. This exactly solvable point is characterized by a pure 
Gaussian decay of the pair correlation. The evaluation of the leading term 
of the (r — 2) expansion for h{r) indicates the change from the monotonic to 
oscillatory behavior just at F = 2 [5j. 

The long-range tail of the Coulomb potential, which is reflected in the sin- 
gular behavior of the Fourier component (11. 5p around k = 0, causes screening 
and thus gives rise to exact constraints (sum rules) for the charge correla- 
tion functions (see review [H]), like the zeroth- and second- moment Stillinger- 
Lovett conditions [71 18]. Their derivation can be based on the exploration of 
the Ornstein-Zernicke (OZ) equation 

h{r, r') = c(r, r') + J dV c(r, r")n{r")h{r" , r') (1.9) 
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relating the pair correlation function h with the direct correlation function c. 
Within the diagrammatic scheme of the renormalized Mayer expansion [1], 
the direct correlation function of the bulk jellium is expressible as 

c(r) = -/5g20(r) + Qeg(r), (1.10) 

where c^eg denotes contributions of all completely renormalized Mayer dia- 
grams. Since these contributions are short-ranged, the Fourier transform of 
Creg has an analytic fc-expansion around k = 0. Consequently, as A; — > 0, 

^ik) = -^ + 0{l). (1.11) 

Writing the OZ equation (11.91) in the 2D Fourier space 

h{k) = c{k) + 2Tmc{k)h{k), (1.12) 

the small-fc expansion of c (11. lip fixes the zeroth and second moments of h{r). 
In terms of the two-body density, these sum rules read 

y"dV [n(^)(r,0) -n^] = -n, (1.13) 
/ dV|rp n(2)(r,0) -n^ = . (1.14) 

It is clear from the derivation procedure that the consideration of a short- 
ranged, e.g. hard core, potential in addition to the Coulomb potential does 
not alter the results f ll.l3p and (11.140 . We add for completeness that for the 
2D jellium also the fourth-moment condition ^ (related to the availability 
of the exact equation of state) and the sixth-moment condition [10] (derived 
within a classification of renormalized Mayer diagrams) are known. 

In this paper, we study a typical situation in the theory of colloidal mix- 
tures [II1[I2]: a "guest" particle with charge Zq is immersed into the bulk 
interior of a Coulomb system, in our case the jellium. Possible values of 
the parameter Z are restricted as follows. When q is the elementary charge 
e of an electron, Z is the valence of an atom and as such it should be an 
integer. In general, the jellium can be composed of multivalent charges 
(g = ±2e, ±3e, . . .) and in that case Z can take rational values. In the 
considered case of the pointlike guest charge and two spatial dimensions, the 
value of Z is bounded from below by a collapse phenomenon. Namely, the 
Boltzmann factor of the guest charge Zq with a jellium charge q at distance 
r, r^^ , is integrable at small 2D distances r if and only if 

Z>-^. (1.15) 

This is the region of the thermodynamic stability for the jellium system plus 
the guest charge Zq. 
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The aim of the present paper is to extend the Stilhnger-Lovett sum rules 
fll.l3p and fll.141) to the presence of the guest charge Zq in the bulk jellium. 
For this purpose, we introduce "conditional" densities: let n(r|Zg, 0) be the 
average density of jellium particles at point r induced by a pointlike charge 
Zq placed at the origin 0. The corresponding induced charge density will be 
denoted by p(j\Zq, 0) = g[n(r|Zg, 0) — n]. Evidently, if Z = 1, i.e. the fixed 
particle has the same charge as the species forming the jellium, it holds 

n(2)(r,0) = n(r|g,0)n(0). (1.16) 

The sum rules f ll.l3p and fll.l4p can be thus rewritten in the form 

y'dVp(r|g,0) = -g, (1.17) 

/dV|r|V(r|g,0) = — (1.18) 
J irTn 

The zeroth-moment condition (11.170 reflects a trivial fact that the charge q 
is screened by a cloud of the opposite charge —q. The condition (11.181) tells 
us that the second-moment of this charge cloud has a prescribed value. Our 
task is to generalize these exact constraints for the conditional charge density 
p{r\Zq,0), where the guest-charge parameter Z lies in the stability region 
(ll.lSp . We notice that there exists one trivial case Z = 0, for which the 
obvious equality n(r|0, 0) = n implies that all charge moments vanish, 

y"dV|r|2-'p(r|0,0) = for j = 0, 1, 2, . . .. (1.19) 

The generalization of the zeroth-moment relation p.l7p is straightforward: 

J dVp(r|Zg, 0) = -Zq, (1.20) 

i.e., the guest charge Zq immersed in the jellium is screened by an excess 
cloud of jellium particles carrying exactly the opposite charge —Zq. 

The generalization of the second-moment relation fll.l8p is nontrivial. We 
would like to emphasize that the derivation of the sum rule fll.l4p . or its equiv- 
alent fll.l8p . using the OZ equation was based on the translation-invariance 
property of the bulk jellium. The introduction of the guest charge Zq with 
Z ^ 1 breaks the translation symmetry of the jellium and one has therefore 
to apply other more sophisticated approaches. Here, we use a mapping tech- 
nique of the 2D jellium with the coupling constant F = (even positive integer) 
onto a discrete ID anticommuting-field (fermion) theory, introduced in ref. 
|13j and developed further in refs. [TU dSl [IS] . The general formalism of the 
mapping technique is briefly recapitulated in Section 2. 

The present application of the fermionic mapping to the thermodynamic 
limit of the jellium in the disc geometry, with the guest charge Zq fixed at the 
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[1.21] 



disc center, is the subject of Section 3. Within the fermion representation, 
a couple of constraints for fermionic correlators is derived by using specific 
transformations of anticommuting variables. Under the assumption of good 
screening properties of the jellium system, these fermionic constraints imply 
the electroneutrality sum rule (ll.20p and the desired second-moment sum 
rule: 

/dV|r|V(r|Z,,0) = --i-Z,[(2-£)4. 

valid in the guest-charge stability region fll.lSp . Although this relation was 
obtained for the series of discrete values of the coupling constant F = 2, 4, . . ., 
we assume its validity for all real values of F corresponding to the fluid regime. 
It is easy to verify that the formula (11.211) is consistent for Z = 1 with the 
result fll.lSp and for Z = with Eq. (11.191) . In contrast to the zeroth-moment 
condition (11.201) . the second- moment sum rule (11.211) provides some exact 
information about possible sign oscillations of the charge cloud screening the 
guest particle Zq and this topic is also discussed in Section 3. 

The exact weak-coupling DH analysis of the studied guest-charge problem 
is presented in Section 4, with the final result 

/dV|rPp(r|^g,0) = asF^O. (1.22) 
J vrFn 

The crucial formula (ll.2ip evidently passes this test. 

The exact treatment of the problem at the free fermion point F = 2, 
performed in Section 5, leads for stable Z > —1 to the result 

/ dV|r| V(r|Zg, 0) = _^^i^±ll at F = 2. (1.23) 
J znn 

The formula (ll.2ip passes also this test. 

Some concluding remarks are given in Section 6. 



2 General formalism 

Let us consider the jellium consisting of N mobile particles confined to a 
2D domain A; the plain hard walls surrounding A do not produce image 
charges. In terms of the complex coordinates {z,z), the potential energy of 
the particle-background system is given by 

E = Eq + q^(l){zj, Zj) - ^ In \zj - Zk\. (2.1) 

j j<k 

Here, 0(2, z) is the one-body potential induced by the background plus per- 
haps some additional fixed charges and Eq is the (background-background. 
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etc.) interaction constant which does not influence the statistical averages 
over particle positions and therefore will be omitted. The canonical partition 
function at the inverse temperature (3 reads 

where the one-body Boltzmann factor w{zj,Zj) = exp[—(3q(f){zj, zj)]. The 
particle density (11. 6p can be obtained in the standard way 

n{z, z) = w{z, ^)^^7^- (2.3) 
dw[z, z) 

For the coupling constant F = 27 (7 = 1, 2, . . . an integer), it has been 
shown in ref. [13] that the partition function (12. 2p can be expressed as the 
integral over two sets of Grassmann variables {C,j"\ ^j"''} each with 7 compo- 
nents (a = 1, . . . , 7), deflned on a discrete chain of N sites j = 0, 1, . . . , — 1 
and satisfying the ordinary anticommuting algebra [17], as follows: 

= JvijVCexp[S{^,i:)], (2.4) 

7(JV-1) 

S{^,^) = (2.5) 

j,k=0 

Here, VifjV^ = Ufjo^ d^j^"'^ . . . dipfd^^"'^ . . . d^f and the action S involves 
pair interactions of "composite" operators 

N-l N-1 

E E (2.6) 

Jl.---J7=0 ki,...,k^=0 
Ul + ---+jj)=j (fcl+--- + fc7)=fc 

The interaction strength is given by 



^zw{z,z)zh^; = 0,l,...,7(Ar- 1). (2.7) 



Using the notation (■■■)= / P^/^P^e"^ ■ ■ ■ jZjq for an averaging over the an- 
ticommuting variables with the action (12. 5p . the particle density (12. 3p is ex- 
pressible in the fermionic format as follows 



7(7V-1) 

'1 ■ 



n(z,z] = w[z,z] E {'^j'^k)z^z''. {2.t 



j,k=0 



Speciflc constraints for the fermionic correlators (Ej'^k) follow from the 
fermionic representation of the partition function as the results of certain 



7 



transformations of anticommuting variables which maintain the composite 
nature of the action fl2.5p . 

Let us first rescale by a constant one of the field components, say 

ef^/i^f^ J = 0,1,...,N-1. (2.9) 

Jacobian of this transformation equals to /i^ and the fermionic action 5* 
transforms to fiS. Consequently, 



= /i-^ / V^jjV^exp ^,^jk^k . (2.10) 

\ j,k=0 ' 



Zj^ is independent of fi and so its derivative with respect to is equal to zero 
for any value of /i. In the special case = 1, the equality c?^lnZAr|^=i = 
implies the constraint 

7(iV-l) 

Y: w,k{^,^k) = N. (2.11) 

j,k=0 

Let us now consider another linear transformation of all ^-field compo- 
nents 

ef^A^ef^ j = 0,l,...,iV-l; a = l,...,7. (2.12) 

Jacobian of this transformation equals to \'y^(^~'^y^ and the fermionic action 
S transforms to Z^J^^g^'' X^EjWjk'^k- Consequently, 

Zr, = \-^''^''-'^/' V^pV^expi ^ A^S,^i;,A . (2.13) 

V j,k=0 J 

The equality d\\rL Z]\f\x=i = implies the following constraint 

7(7V-1) , 

Y j^,fc(S,vl/,) = -7iV(iV-l). (2.14) 

j,k=0 ^ 

The application of the transformation fl2.12p to all ip-Held components leads 
to the complementary condition 

7(iV-l) . 

J2 kw,k{^,^k) = -lN{N-l). (2.15) 

j,k=0 ^ 
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3 Derivation of sum rules 

We study the jellium model confined to the domain of disc geometry A = 
{r, r < R}, with the guest charge Zq fixed at the origin 0. The guest charge 
Zq together with the total charge Nq of mobile particles are compensated 
by the fixed background of charge density —nq via the overall neutrality 
condition 

Z + N = TTR^n. (3.1) 

The potential induced by the homogeneous background is g7rnr^/2, the guest 
charge interacts with jellium particles logarithmically —Zqlnr. The total 
one-body potential acting on each particle 

2 

0(r) = g^^^ - Inr (3.2) 

possesses the circular symmetry. 

At the coupling F = 27 (7 = 1, 2, . . .), the one-body Boltzmann factor 
w(j) = exp[— /30(r)] reads 

w{r) = r^'*'^ exp(— 77™r^). (3.3) 

Within the fermionic representation of the jellium fl2.4l) - fl2.8p . the interaction 
matrix (12. 7p becomes diagonal 

Wjk = SjkWj, Wj = / d^r r^^'''^"*"-'-' exp(— 77rnr^). (3.4) 

The consequent diagonalization of the action (12.51) in composite operators, 
S = J2]^~^^ '^j'^i^ implies that = 5jk(^j^ j) and the representation 

of the particle density (12.81) simphfies to 

7(iV-l) 

n(r|Zg,0) =e-^™' ^ (S,-^,)r2(^^+^). (3.5) 

i=o 

The constraint (12.111) is expressible as 

7(iV-l) 

^ w,{E,^,) = N (3.6) 
and the couple of complementary conditions (12.141) and (I2.15P reduces to 

7(7V-1) . 

^ 3w,{'^,^,) = -^N{N-l). (3.7) 
3=0 
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Using the definition of tlie interaction integrals fl3.4p . it is easy to sliow 
that the constraint fl3.6l) is equivalent to the relation 



d^r n{r\Zq, 0) = A^, 



(3.8) 



which reflects a trivial fact: the total number of mobile particles in the disc 
domain A is equal to N. With regard to the electroneutrality condition (13.11) . 
the relation (13.81) can be rewritten in the form 



By a simple analysis we shall argue that this condition involves in fact two 
sum rules, the bulk one and the surface one. Let us divide the disc do- 
main A onto its "bulk" part A^ = {r, r < R/2} and the "surface" part 
As = {r,r = R — X with < a; < R/2} [x denotes the distance from the disc 
boundary) and rewrite Eq. (13. 9p as follows 

l-R/2 rR/2 



27rrdrp{r\Zq, 0) + / 2tt{R - x)dxp{x\Zq, 0) = -Zq. (3.10) 



Let us assume that the system of charges has good screening properties, i.e. 
the decay of particle correlations at large distances r is faster than any inverse 
power law, say exponential oc exp(— Kr) with k being the inverse correlation 
length (like it is in the weak-coupling limit F ^ 0) or even Gaussian oc 
exp[— (/tr)^] (like it is at the free-fermion point F = 2). In the R ^ (yo limit, 
the particle density differs from the constant n only: in the bulk region close 
to the disc center (up to r ~ and in the surface region close to the 
a; = boundary (up to x ~ k^^). The charge profile close to the boundary 
p{x\Zq,0) is influenced by the screened guest charge Zq (exponentially or 
even Gaussianly) weakly as i? ^ oo. Forgetting these small terms, one can 
put 



pix\Zq, 0) ~ p{x\0, 0) = pix) + -h{x) + — /2(x) + ■ ■ ■ , (3.11) 



where the long-ranged inverse-power-law terms 1/-R, 1/-R^, . . . are due to the 
nonzero curvature of the disc surface and the respective coefficients /i, /2, • • • 
are short-ranged functions of the dimensionless parameter nx. Thus, Eq. 
(I3.10p splits in the limit R ^ oo into the Z-dependent bulk electroneutrality 
condition of present interest 




(3.9) 





1 




(3.12) 



and a series of Z-independent surface conditions 




(3.13) 
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the lowest one of which takes the form of the surface electroneutrahty 

dxp{x) = 0. (3.14) 



To make use of the constraint (13. 7p . we first differentiate both sides of the 
density representation (13.51) with respect to r, then multiply the result by r 
and finally integrate over the disc domain, to obtain 



[ d'^rr—n{r\Zq,0) = 2'jZN-2'j7m [ d^r r'^n{r\Zq,0)+2 V jWj{Ej<iJj). 
or JA 

(3.15) 

The Ihs of this relation can be integrated by parts, the summation on the rhs 
is given by the constraint of interest (13.71) . After simple algebra, the relation 
(I3.15P is transformed to 

— 2TC'yn / d^r r^p(r|Zg, 0) = {2 — 'y)Zq + 'yZ'^q 

J A 



+2t^W 



p{R\Zq,Q) + ^qn 



. (3.16) 



Like in the previous analysis of Eq. (13.91) . we divide the disc domain A onto 
its bulk and surface parts to express the integral in Eq. (13.161) as follows: 



/ 27rrMrp(r|Zg,0) + / 2'k{R - xfdxp{x\Zq,0). (3.17) 
Jo Jo 

Under the assumption of good screening properties of the jellium, the bulk 
and surface regions are coupled weakly in the R ^ oo limit and one can 
consider once more the expansion (13. lip for the boundary charge density. In 
this way, one gets from Eq. (13.160 the Z-dependent bulk condition 



rR/2 



[ dV r'^pMZq, 0) = —Zq [(2 - 7) + jZ] 

J 2TT'jn 



(3.18) 



which is equivalent after the substitution 7 = T /2 to the one of primary 
importance (ll.2ip . and a series of Z-independent surface conditions 



2Tt-fn / 27r(i? - a;)^dxp(x|0, 0) = 2nR^ 
Jo 



p{x = 0|0,0) + ^qn 



. (3.19) 



The lowest-order surface condition can be obtained by summing Eq. (13.130 . 
multiplied by 2TT'ynR'^, with Eq. (13.190 . The final result reads 



p{x = 0) = "■^Q'^ + 477771 / dxxp{x). 



(3.20) 
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This relation is known as the contact theorem [TSl [T21 ^U\. Although all 
relations were derived for 7 = r/2 a positive integer, it is reasonable to 
extend their validity to all values of F corresponding to the fluid regime. 

As was mentioned in the Introduction, the generalized second-moment 
sum rule (11.211) is consistent with the available results (11.181) for Z = 1 and 
(I1.19P for the trivial case Z = 0. In the next two sections, we test this sum 
rule also in the weak-coupling F ^ limit (Section 4) and at the free-fermion 
point F = 2 (Section 5). 

In contrast to the zeroth- moment electroneutrality condition (I3.12p . the 
generalized second-moment sum rule (I3.18p . or equivalently (ll.2ip . provides 
an exact information about possible sign oscillations of the induced charge 
density p{r\Zq,0) in space. If Z > 0, the guest particle and jelhum charges 
repeal each other and therefore p{r\Zq,0) ~ —qn as r — > 0. Provided that 
p{r\Zq,0) does not change the sign when changing r from to 00 (where p 
vanishes), its second moment has the sign opposite to Zq. Similarly, if Z < 0, 
there is an attraction between the guest particle and jellium charges, so that 
p{r\Zq, 0) goes to infinity as r — * 0. Consequently, when p{r\Zq, 0) does not 
change the sign when going from r = to r ^ cxd, its second moment has 
again the sign opposite to Zq. The sufficient condition for sign oscillations of 
the charge density p{r\Zq, 0) in space is that its second- moment has the sign 
of Zq. In view of the result (ll.2ip . the sufficient condition for oscillations is 
that the guest-charge parameter Z lies in the interval 

- ^ < Z < 1 - 1, (3.21) 

where the lower bound, see Eq. (11.150 . ensures the thermodynamic stability 
of the pointlike guest charge Zq. The inequalities (13.210 have no solution for 
F < 2. For F > 4, there exists also an interval of positive values of Z for 
which the induced charge density certainly exhibits sign oscillations. 



4 Weak-coupling limit 

The effective potential at distance r from the guest charge Zq, placed at the 
origin and surrounded by mobile g-charges of the average density n{r\Zq, 0) 
plus the neutralizing background of charge density —qn, is given by the 2D 
Poisson equation 

A0(r) = -27rg {Z5{r) + [n{r\Zq, 0) - n]} . (4.1) 

The weak-coupling (high-temperature) region F — >^ is described rigor- 
ously by the Debye-Hiickel theory [31 |2T] . Within this mean-field approach, 
the average particle density at a given point is approximated by replacing the 
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potential of mean force by the average electrostatic potential at that point, 
n(r|Zg, 0) = nexp[—(3q(f){r)]. The mean-field Boltzmann factor can be lin- 
earized at high temperatures, exp[— /3g</)(r)] ~ 1 — jSqcpir). The Poisson Eq. 
(14.11) then reads 

(a - K^) 0(r) = -27rZg5(r), (4.2) 

where k = \/2'KTn is the inverse Debye length. 

Due to the circular symmetry of the problem, A = d'^ + {l/r)dr. Eq. 
(14.21) . subject to the condition of regularity at r — oo, thus implies 

0(r) = ZqKoM, (4.3) 

where Kq is a modified Bessel function |22j . 

The induced charge density around the guest charge Zq is obtained in the 
form 

p{r\Zq, 0) = -ZqnTKo{Kr). (4.4) 

Since the stability lower bound (ll.lSp is Z > — oo in the limit F — > 0, this 
result applies to all real values of Z. The charge density (14.41) is always 
a monotonic function of the distance r which keeps its plus {Zq < 0) or 
minus {Zq > 0) sign in the whole interval of r G (0, oo). Its moments Ij = 
/q°° p{r\Zq, 0) {j = 0, 1, . . .) are given by 

roc / 2 \ ^-^ 

Ij = -ZqK^ j drr'^^+^Ko{Kr) = -Zq i^-j [r(l + j)] , (4.5) 

where T{x) denotes the Gamma function. For j = 0, the electroneutrality 
condition (11.201) takes place. For j = 1, one arrives at the second-moment 
formula (11.221) which is in full agreement with the general result (11.211) taken 
in the weak-coupling limit F ^ 0. 



5 The free-fermion point 

The fermionic representation of the 2D jellium simplifies substantially for 
the coupling constant F = 2 (7 = 1), because the composite variables (12. 6p 
become the ordinary anticommuting ones. Having the fermionic action of the 
form 5* = Y.f=o^ ^j^j^j easy to show that 

Af-l 

Zm = Yi ^i' (5-1) 

j=0 

{^j^j) = — j = 0,l,...,iV-l. (5.2) 

Wj 
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In the limit of the infinite disc radius i? — oo, the interaction strength fl3.4p 
at 7 = 1 is given by 



1 



^^ = Z7ZZ^^nZ + 3 + l). (5.3) 



For an infinite number of jelhum particles N 00, the particle density (13. 5p 
induced by the guest charge Zq reads 

n{r\Zq,0) ^ , 2^ , u^ -t^ t^^' 

-^^±-lL-L = fUnrn, /^(t) = e*> — -. (5.4) 

n ^ ^ ' PoT{Z + j + l) ^ ' 

It is seen that the induced density is well defined for Z > — 1, and this is 
indeed the range of the guest-charge stability (11.151) for F = 2. 

Let us first treat the region of Z > (g > will be considered for 
simplicity). We shall need the incomplete Gamma function which is defined 
as follows |22j : 

/•oo rt 

T{Z,t)= dss^-V = F(Z) - / dss^-V^ Z>0. (5.5) 



It can be readily shown by applying a series of integrations by parts that 

00 .Z+j 

r(z,t) = r(z) -r(z)e~'y — -. (5.6) 

The function fz(t), defined in Eq. (15.41) . is therefore expressible as 

Mt) ^ 1 - (5.7) 

and the induced charge density reads 

/ , N F(Z, vmr^) 

pir\Zq, 0) = -qn^,^^^, Z > 0. (5.8) 

Since dtT{Z,t) = —t^~^e~^, the derivative drp{r\Zq,0) is positive for any 
value of r. Consequently, the induced charge density is the monotonically in- 
creasing function of r, going from —qn at r = to at r ^ oo. The moments 
of the charge cloud around the guest particle Ij = 2TTrdr r"^^ p{r\Zq,0) 
(j = 0, 1, . . .) are given by 

2,-+i F(Z,W) q T{Z + j + l) 

Ij = -2nqn drr ^+ — — — — = - — — — , 5.9 

Jo T{Z) {j + l){7my T{Z) 

where we have applied an integration by parts. For j = 0, one recovers the 
electroneutrality sum rule (ll.20p . For j = 1, one gets the result (11.230 which 
is in full agreement with the general result (11.211) taken at F = 2. 
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As concerns the stability region of negative Z- values — 1 < Z < 0, we first 
write down a recursion relation for fz{t) following from the definition fl5.4p : 



fz{t) = e 



-t 



+ fz+iit). 



(5.10) 



T{Z + l) 



Thus 



p{r\Zq, 0) = qn e 



—nnr- 



.2 (vmr^)^ 



T{Z+l,7rnr^) 

r{z + i) 



Z>-1. (5.11) 



T{Z + 1) 



The induced charge density is now the monotonically decreasing function of 
r, going from oo at r = to at r ^ oo. It is easy to verify that the formula 
for its even moments coincides with the previous one (15. 9p . This fact permits 
one to extend the validity of the zeroth-moment (11.201) and second-moment 
(I1.23P sum rules to the region of negative Z-values —1<Z<0. 

6 Conclusion 

In this paper, we have generalized the standard zeroth- and second-moment 
Stillinger-Lovett sum rules for the charge correlation functions to the pres- 
ence of a guest charge immersed in the bulk interior of the 2D jellium. The 
derivation procedure was based on the fermionic technique which is associ- 
ated specifically with the 2D jellium model. It is an open question whether 
the generalization of the sum rules can be accomplished also in higher di- 
mensions or for many- component Coulomb fluids. The present results might 
inspire specialists to establish some new phenomenological arguments which 
go beyond the standard ones. 
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